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FURTHER COMMENTS ON THE APPLICATION OF THE METHOD 
OF AVERAGING TO THE STUDY OF THE ROTATIONAL MOTIONS 
OF A TRIAXIAL RIGID BODY, PART 3 

1. Introduction 

In [J.R.,1970] , we derived a set of differential equa- 
tions, analogous to the Lagrange planetary equations, for use 

0 

in studying the perturbations of the rotational motions of ex- 
tended rigid bodies. These differential (variational) equations 
are expressed in terms of the perturbing torques and they are 
valid for either conservative or nonconservative torques. In 
[A. R. ,1971], [F.R.,1972] and [M.R.,1973], these variational 
equations were applied to the case of a rapidly spinning tri- 
axial body moving in an elliptic orbit, in which the orbital 
plane is regressing at a constant rate. The explicit differ- 
ential equations obtained in this application were integrated 
by the method of averaging (described in [F.R.,197l]) to develop 
secular analytical express ions, which, to first-order In a 
small parameter, describe the complete space motions of the 
rigid body under the influence of nonresonant gravity-gradient 
perturbations. 

In this report, the effects of aerodynamic torque on the 
rotational motion of an orbiting satellite will be studied, as 
another example of the application of the variational equations 
derived in [J.R.,1970] and the method of averaging discussed 
in [F.R.,197l]. As an artificial satellite moves through the 


References to our earlier reports of June 16, 1970, 
February 19, 1971, August 2, 1971, February 21, 1972, and 
March 27, 1973 are indicated by [J.R.,1970], [F.R.,197l], 
[A.R.,197l] , [F.R.,1972]and [M.R.,1973]> respectively. 
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Earth' s atmosphere* there is a drag force introduced by the 
interactions of the oncoming molecular flow with the satellite 
hull. This drag force* which is one of the most important of 
the perturbing forces influencing the rotational motion of the 
satellite about its center of mass* brings about secular var- 
iations in all of the variables, i.e., ^„* 9„, <p ^Q'* <£'*h. 

n. H H T 

In particular* and unlike the gravity-gradient force* this drag 
force brings about secular variations in both the magnitude of 
the angular momentum vector h and the coning angle 6' . 

Beletskii [1 ] assumes a plastic or perfectly inelastic 
impact of the molecules of the upper atmosphere with the sate- 
llite hull and neglects the effects due to the rotation of the 
Earth's atmosphere. He further assumes that the linear speed 
of the satellite hull with respect to the center of mass is 
very small compared to the speed of the center of mass itself. 
He then derives a symbolic expression for the aerodynamic tor- 
que N a about the center of mass* for a triaxial body with a 
general exterior surface. His expression is valid through terms 
of first order in the magnitude of the angular velocity of the 
rotation of the satellite relative to the atmosphere. His as- 
sumptions will be used in this report in deriving an explicit 
expression for the aerodynamic torque in terms of a convenient 
cylindrical coordinate system. This general expression will 
reduce to the expression for a satellite with a surface of 
revolution given in [1]. 
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The explicit expression for the torque for a satellite 

with an arbitrary shape involves surface integrals whose limits 

depend upon the instantaneous surface of attack. Although this 

torque expression is analytically intractable for the general 

case* it does become manageable if the satellite has a surface 

of revolution. For this reason, the considerations here are 

restricted to the study of the secular effects of aerodynamic 

torque on an orbiting, uniaxial satellite. Beletskii [1] uses 

-y- 

the same torque expression to study the secular motion of the 
angular momentum vector for a rapidly spinning, uniaxial, orbit- 
ing body with a surface of revolution by applying the method of 
averaging to the Euler's dynamical equations. In this report, 

however, the secular variations of all the six variables which 
describe the complete rotational motion of such a body about 

its center of mass under the influence of aerodynamic torque 
are studied by applying the method of averaging to the varia- 
tional equations. Approximate, averaged, first-order, ordinary 
differential equations for the six variables are obtained, 
which are particularly convenient for numerical computation. 

If the problem is further simplified to the case of an uni- 
axial body moving in a circular orbit with constant air density, 
the approximate, averaged differential equations for the magni- 
tude of angular momentum vector and coning angle can be Inte- 
grated. It is found, as a verification of the results given 


There is a misprint in the torque expression for a sate- 
llite with a surface of revolution given in [1] [equation (1.3. 11 
p. 16] . In a private communication, Beletskii agrees that the 
correct expression for the torque is given by equation ( 3 * 17 ) > 
of this report. 
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in [1], that the magnitude of the angular momentum vector de- 
creases exponentially with time approaching zero as a limiting 
value# while the coning angle may either decrease or increase 
with time approaching either zero or ot/2 as limiting values, 
respectively. 

2. Coordinate Systems 

With reference to Figure 2.1, let Ox* y‘ z‘ represent a 
body-fixed system whose origin 0 coincides with the center of 
^hn'ass of the satellite. The rectangular coordinate system 

. _ .a, 

0x‘ ‘ y‘ 1 z‘ 1 is obtained by rotating the system Ox’y'z 1 about z* 
through the angle f3 v in the sense shown. The angle is, 

determined by the condition that the y'*-axis lies in the plane 
determined by the z'-axis and the velocity V q (translated to 0) 
of the center of mass with respect to an inertial frame at geo- 
center. The x 1 ‘-axis is chosen to form a right-handed system. 
Let i.' , j', k‘ and i_‘ ' , j‘ * , k‘ * , be the unit vectors along the 
x*, y* , z 1 and x‘ ’ , y‘ 1 , z‘ * axes, respectively. The angle & v 

is drawn from the z‘ -axis to the line determined by V . The 

— o 

cylindrical coordinates z* of an arbitrary surface point P 

of the body is introduced where ~y is reckoned from the y‘ -axis 
in a counterclockwise direction in the x'y' -plane. The unit 
vectors associated, with the cylindrical coordinates ^ , y , z‘ 
are designated by e , e , e , , 

P r 2 


respectively 
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The following relations can be obtained readily 


e v = s c J' 1 + c f k' 1 , 
v -iyf— d v 



( a) 





Vv 

0 

1" 

= 

1 

CO 

< 

Vv 

0 

k* ' 


0 

0 

1 



i’ 


(b) (2.1) 


k* 



0 

0 

1 


i* 

y 

k' 


( c) 


where e^. represents the unit vector along V . 


3* Aerodynamic Torque 

Assume that the impacts of the air molecules of the 
oncoming flow with the surface of the satellite are plastic or 
perfect inelastic and neglect the effects due to the rotation 
of atmosphere. Beletskii [1] has shown that the aerodynamic 
torque N , with respect to the body-fixed system, through 

A 

terms of first order in the magnitude of the angular velocity 
of rotation of the satellite relative to the atmosphere, is 
given by 
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5a = I °fa V ° ^ (s --V )(a v X£) 


dS 


(3.1) 


+ 1 c <0 V 


fi V C . [.(n.cox r) ( e,r x r)-+(n*eJ-(oJxr) xrldS 
2 (a ° ) - -V V - - 


where c is a coefficient. The quantities p> n, r, and uJ 
represent the atmospheric density, unit outward normal of the 
surface element dS of the satellite, the position vector of 
the surface element dS and the angular velocity vector of the 
satellite's rotation about its center of mass, respectively. 
{The domain of the surface integral is indicated by S and it 
is defined by the Inequality 


I Q • n i' ° (3.2) 

The remainder of this section will be devoted to the 
development of an explicit expression for the torque (3.1) in 
terms of the cylindrical coordinates p , y , and z 1 . The unit 
outward normal to the surface element dS can be written either 
in the component form 


or 


n = n ^ e ^ +1 n e + n , e « 

“ |° — 7 ~y z 1 £z' 

= n x , V + n yl £ + n z , k' 


( 3 . 3 ) 


n 


( 34 ) 
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Using equation (2.1(c))* we can write 


n , = - P o~ s -1 n c v 

x ' r n y — ar ar 


n 


P r n V - V-V'r 


( a) 


(b) 


(3.0 

r 


where n n = p (T • 

( [n 

Let cf>( , y > z*) be the relation which describes the 
surface of the satellite. It is found that 


e + 


W = <* f e f ^ 


4> , e (3.6) 

’ 7. ' 7. * 


and that 


n = Z F <Pp*e fi +a/ r W f e r +y 2 ,a z . (3>?) 

“ * * P 3 : 172 

fop + <?• Z ,1 


^ ^ ,<$)■= 2 0 
i° ? < f V /° ^ 


(3.8) 


<Py = gg > f Z l = ^ 

9 i ^ z 1 

r 

Comparing expressions (3*3) and (3*7) * we have 


n = 

P 


2p £ 


9' 


2 2 


-2 2 


H * 2 V2 

fo P <? p* + i p > <*> y + z .i 


(a) 



8 


cj> 


y rl - ~s z 2 2 

9 ^2 + { ( o) z .3 


W 


(b) 


Y 


(3.9) 


n 


7 Z . 


r , - 2 _2 2 i 1/2 

+ ( f } + ? Z‘ ] 


( c) 


When referenced to the body-fixed system, the position vector 
of the surface element dS has the following form 


+ y 1 j. 1 + z* k* 


r = x‘ 


where 


7 

/ 


x* = - p a 


r 


r = 


(3.10) 


(3.H) 


Using equations ( 3*5) » ( 3 *9) , ( 3 *10) , ( 3 *11) and the equations 
of transformation (2.1), we can express (3*1) l n the more 
explicit form 

H A = 1 cp V “( N A1 , i' + N A1 , J' + n a1z , k') 
d 1 a 


+ 2 C P V (J 11^ + J 12^y' + ^ 


+ < J 21^x‘ + J 22 a) y‘ + J 23 a3 z'^ 


(3.12) 


+ ( J oo + J oO + J oJ ) k 1 ] 
31 x' 32 y* 33 z' “ 
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where, u) x ,,o) and go t are the rectangular components of 

y ^ 

cJ expressed in the Ox'y'z' system, and 


N 


Alx* 


= )(z' s r c. -Pc. c ) n dS 


( a) 


W = ' (z ’ s <- 3 


r s ^ -pc s 
" V V I $ V t 


c s ) n dS 
e 


(b) 


N A1z* 


= (ps c c a a , - p s r s „ c ) n dS 

J* I *v Pv x f *v Pv y e 


( c) 


J. 


= 1 i r(n = 


n - j v'-«* ' z ' r ^r, + y \ % 1(21 a i V v v 

s * r 


p °* 


c - c _ ) -n ( Z 1 2 + // c") <r dS 


V ^ 


e 


r v 


(d) 


'12 •' 



[(n . - z 1 f ) p s -1 z ' n c y ](z' s, c . 

2 n i r j- rx P 


V 


-pc, o ) - ik ^ s c ( dS 

’ s v r E I r x ) 


( e) ( 3 . 13 ) 


J = f fn ( z’ s c c - p c e c ) -n z'p l dS ( f ) 
13 J y *v Pv ' r e r r 


21 


[<n z , - •'<r B )f,a r + 1 ]■(«'•, T y 


V 


-pc c s ) -n o ^ s c ( dS 

r Sv Jr ®f )r irj 


(s) 
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22 


= [ { f(n z , - z><T n )p 


« - z' fjps -I. 2 1 n c y ]{z* s r s 




f 


r r 


5 v ^ v 
(h) 


■ r° r vV ) ’ ne( *‘* + ds 
r a 3 = / fn y u ‘ s f v > v -p-°f v V- + VY °r ]dS (1) 


J 31 = 



[(n 2 , - zT n )(>c y «• n. e 


v a v y 


T'rrfl''** ‘YVI as 


( j) 


/32 




|t(n z , - - i- *’ » r V H f S S 


c . s 

v Z 3 v ^ 


'? °» j 


psc s „ c ) + n z' p c ( dS 

( s v Pv r 9 


= 5 [n( f 


(k) 


(3.13) 


c C s - /? a s c) -n /> ] dS ( jt ) 

33 ^ , I *v-hi I Ivfvi' e r 

s< 


n e = » <■ 3 „ ( i? <r 3 + 1_ n e ) 

^ v r v ’ n 0 p 2T ? 


+ s s c ( £> (T c r -1 n s ) + n , c r 

s vPv' nf t r y z *v 


(m) 


Equation (3*12) gives the aerodynamic torque N A for a 
body with a general surface. Even for short periods of time 
it is observed that N is time dependent> principally through 

n 

and ^ y. For longer periods of time V Q and u) must also be 
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considered to be time dependent vector quantities as they ap- 
pear in N^. The expression becomes simpler for a body with a 
surface of revolution with respect to z'-axis since n y= 0* 
Moreover# if the semi- cons trained system Ox’ 'y'^ 11 is chosen 
as the reference system# the simplications 


P V = °* 03 X’ =cu x" 


f°* dr = f •» c >- d ^ = 


r' 2 

s c d^- = 0 


* y 


(3.14) 


may be Introduced and then it is found that 


N A1x« = N Aly« " J 1 2 = J 13 = J 21 ~ J 31 = ° 


(a) 


N A1z' = I 5 ( s J o v x k 


V -V 


(b) (3.15 


J 11 = " I 5 " -c J y J ^ z ' 2+ j> )n z* + (n z ,_ ^n 21 ' V 


/“ 2 2 2 
-8 f f( Z 1 + p C 


f c f l Vv ,,( v - f /Vv |as 

( c) 


J 22 = 'V " C Sv I [(n2, ' Z ’ f* 3 * + n z' U ’ a+ r 2 v )lds 


' S £ v 5 + f^r ) T nf 0 J- aS 


(d) 
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J 23 = *4 ~ 3 S v J 


a 2 


_X 


‘ ,<r nf dS + 1 ".'f 


J 32 “ I 2 = ° 


s y f* , WV s + *f v f' V p V 


+ f ^ V (n z« " zt<f r) ] dS 


(f) (3.1?) 


j. _ w r f n , p 2 dS + s 1 ^ p^ c y- dS 

33 i £ v 2 r « v j n r r 


= -I, = C 


where 


5( J v )=c(W 1 y + w 2 3 S' y - W 3 ° S v ‘“‘Jv (S) 


W 1 “ 


2 2 


(*' n z' * CTj! f 5 °J ) dS 


(b) 


(3.16) 


W 2 = \ *' ^ nt °t dS 


( c) 


W = ( n /) c , dS 


f V 


(a) 


Then expression (3.12) becomes 


N 


A 


= 1 f a 7 o = ( ? V > % * + I °fa V o 1 


(3.17) 
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where 


i=- V^x” 1" + (- I 3 t O y ,,-H 4 ^ ii ) J" + (I 2 ^ ylI -I 1 a> z ,, ) k> 

(3.18) 


Expression. ( 3 . 17 ) is exactly the same_as given in [1] . . 

The aerodynamic torque (3*17) can be- analyzed in two sepa- 
rate parts. Thefirst part, which arises because the center of 
pressure does not coincide with the center of mass of the sate- 
llite, will be referred as the restoring torque N . The second 

AR 

part which introduces damping effects caused by the spin of the 
satellite with respect to its center of mass, will be reffered 
as the dissipative torque N AD * In the notation just described, 
( 3 .I 7 ) can be rewritten in the form 


N = N + N 
“A “AR “AD 


( 3 . 19 ) 


where 


4r = i fa V o S V> % X k " 


(a) 


( 3 . 20 ) 


-AD 


1 c P 

2 I 


a o — 


(b) 


ij.. Approximations of the Aerodynamic Torque 

Equation (3.12) gives a description of the aerodynamic 
torque for a body of arbitrary shape. In principle, the torque 
is to be calculated by integrating over the surface of attack 
S . It is noted, however, that the integrals which are involved 
may be analytically quite intractable since the integration 
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limits depend upon the surface of attack S . In turn, the sur- 
face of attack may depend discontinuously upon the time through 
the variables ft and $ ^ , simplification may be achieved in 
certain important special cases, however. For example, in this 
report, the perturbation caused by aerodynamic torque will be 
analyzed in the special case where the satellite possesses a 
surface of revolution and the semi- constrained system is chosen 
to be the reference. In what follows, the considerations are 
restricted to this special case. 

It is indicated by Beletskiifl] that the principal quanti- 
tative and qualitative effects of the "^aerodynamic torque which 
are common for various bodies of the type considered, can be 
described by using certain approximate formulas for and 1^. , 
where i= 1,2 ,31 j=l ,2, 3,4*5* These formulas are described 
briefly in this section. 

( a) Restoring Torque 

First of all, attention is concentrated on the restor- 
ing torque which is free of the spin of the satellite. Construct 
surface S^, which is parallel to unit vector and S Q , which 
is perpendicular to e^, in such a way that if combined with 
the surface of attack S , a closed surface in space is formed, 
as shown in Figure 4*1 • 

From equations (3*1) and (3.19), we have 




V e„ x 
o — V 


( n • e v ) r dS 


S 


* 


( 4 * 1 ) 
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Let S + S, + S be a piecewise smooth orientable surface and 
1 o 

let u be a vector function of the coordinates (x'lj 1 , z') which 
is continuous and has continuous first partial derivatives with 
respect to the coordinates in some domain ^containing 
S + S-^ + S . Using the divergence theorem, we may write 


L 


u*n dS = 


S +S- +S 
1 o 


j 7- u 

r 


d r 


( 4 - 2 ) 


and, since 


^ x ' — v * — ^ A* + (y* ®v • e) i' + ( z ' ®v*^ k ' (4-3) 


we can write 


f (u • n) r dS s i» f 7 *(x' e v )dr + J.' f 7 - ( y* e y )df 

K* U r 


s* +S..+S 
1 o 


k' I Y *( z ’ ®y) d7 *- 

V 


Since _e^ is independent of t^e coordinateo x' , y', and z*, it 
can be shown that 

(n ’ e v ) r dS = T (4-4) 


1 


s +s x + s 0 


where T is the volume enclosed by the surface S + 



) 
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Since the surfaces and S Q are so constructed that 


*1 * e v = 0, 


£0 * e v = “ 1 


( 4 - 5 ) 


the expression for the restoring aerodynamical torque takes 
the form 

' dS (4.6) 

where r 1 is the position vector of an arbitrary surface element 
of S Q with respect to the center of mass of the satellite. 

Let O' be the projection of the center of mass on S 0 , r^ 
be the position vector of 0* , and r* be the position vector 
of any material point of S Q with respect to 0* . We have, along 
. surface S Q , 

= L 0 + £g (4-7) 


-AR 


= 1 c 


J. C P 
2 \ 


v o£ V 


p 

V 


and hence 


N ad = 1 c p vf S_( e,, x r' ) 
— AR p f a 0 0 ” v -os 


( 4 - 8 ) 


where r rto is the position vector of the centroid of surface S„ 
— o s o 

with respect to O’ . The vector product e, r x r* will be either 

—V —os 

in the positive or negative sense of x' 1 . The magnitude of 


-AR is 


w 


AR 


= 1 c f> V* S n r* , r* = jr’ I 
( a o o os os |— os I 


( 4 - 9 ) 


As shown in Figure 4*1> r = Z^ s r , where is the 

^ OS O £ y o 
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distance between the center of mass and the center of pressure* 
Z Q = OP. Let the projection of P on S Q be P'* the centroid of 
S Q . Equation (4*8) can be rewritten in the form 


s'lc p V“ (e, r x k') 

“AR p f a ° 9 ° v ~ . 


( 4 . 10 ) 


Comparing (4*10) with expression (3»20(a)), it is found that 


C( S y) - C S 0 ( 


(4-11) 


For perfect inelastic collisions* it is also true that 


®( ^y) = ' ©(it - 5 y) 


(4-12) 


Note* from equations (3»l^(c)) and 3.16(d)), that* at an 
instant of time at which either = 0 or S ^ = n* W^ = W^ = 0, 

This suggests the approximations 


W 2 = 3 f y 


W 3 = ^ v f 3 ( f V ! 


( a) 


(b) 


(4-13) 


Where f ^ and f^ are functions of ^y« 
becomes 


Then equation (3»l6(a)) 


c ( £ y) — C [W, C p + f ^ “ ( f p+f _,) C j ] 


1 * V *2 


2 3 ' $ 


( 4«i4) 


which is a polynomial in c r . We may write (4*14) l n the form 


C( £y) = C 


o + C 2 c <f 


(4*15) 
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where 


°o - °< W 1 Y * i 2 ) 

(a) 

(4.16) 

C 2 = ~ c (f*2 + ^'^) 

(b) 

In general C Q and C ^ depend slowly upon time. 

A meaningful 


simple approximation to (4«l4) is obtained by assuming that 

both C and C_ are constants which can be determined either by 
o 2 

using equation (4»H) or any alternative method such that 
expressions ( 4 « 15 ) give reasonable approximation of c( S’ ) • 

If the collisions between the molecules of the oncoming 
flow and the surface of the satellite are not perfectly in- 
elastic and reflections occur# then equation (4»12) does not 
apply. In this circums tance, the simple approximation 

g 

c( S y) = + C x C ^ + C 2 c ^ ( 4*1?^ 

where = c = c W^, may be introduced. If however# 

the satellite is also symmetric about a plane which is perpen- 
dicular to the sxis of symmetry, equation (4«12) will hold ir- 
respective of the reflections of the molecules and equation 
(4.15) will describe the case approximately. For a small angle 
of attack, we may even choose 

c = Cq + C + C = constant (4»18) 

( b) Dissipative Torque 

Next consider the problem of approximating the dissipa 
tive part of the aerodynamic torque given by equation (3.20(b)) 



Approximate formulas will be assumed for 1^, j=l,2,3 >b*5* 'Re- 

.calling the reasoning for equations ( if. 13) * it may he observed 

from equation ( 3 . 1 5 ) that, whenever the factor c , appears 

of 

in an integrand, the integral may be assumed in the form 
8 £ ^ v } ‘ Therefore, we may consider that 1^, 1^ and 

are positive quantities for any value of <P and that the 
principal parts of these functions are constants. It is also 
noted that when either £ ^ = 0 or %, the difference I - 1^ 
is zero. Thus, the principal parts of .1 and 1^ may be approxi- 
mated by the same constant. As to Ig and 1^, they are functions 

of the form s c f .( $ ) . In summary, a simple approximate, for 
* V J v 

may be obtained by assuming that 


*3 " : S = C H 


( a) 


*1 ■ °33 


Ip = C?ps 


32 a ff 


V 


(b) 

( c) 




% = V ! V 


(a) 


where' the C«-, are; constants. Reasonable estimates of these 


constants can be obtained by averaging the values of I. at 
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5» First-Order, Secular Solutions for a Uniaxial Body with a 
Surface of Revolution Subjected to Restoring Aerodynamic 
Torque 

Let the center of mass of the satellite move along an 
elliptic orbit about an attracting center at 0. Assume that the 
inclination angle 0 ° of the orbital plane, the semimajor axis 

o 

"a , the eccentricity e and the rate of orbital precession _Q_ 
are constants and let w represent the true anomaly. We can 
then write the well-known equation from orbital theory that 

V Q = (/t<-/p *) l/2 f, f = (1 + e 2 + 2 e c w ) l/2 (£. 1 ) 

where yU, is a constant equal to the product of the sum of the 
attracting mass and the satellite's mass with the gravitational 

y — 

constant and p = a(l - e 2 ) . Let 0 J ^ ^ 1 represent a 

rectangular coordinate system in which the 5 ' axis is perpen- 
dicular to the orbital plane, positive in the sense shown in 
Figure %•!. "The positive 5 ' axis is chosen to coincide with 
the half line from 0 through perigee. Then if 1 is chosen to 
complete a right handed coordinate system. The unit vectors 
associated with the 0 -f 1 ^ ' system are designated by 

i-*j_.»k_, respectively. 

$' n' “S' 

It is clear that 

i-. = ° w i° * \ J° (a) 

0 

i_, = 8 w i° + <\,i° < b > 

® v = - I ts w i-, -(e + o w ) i_,l (c) 

I 5 ' 1 


(5.2) 
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If the definitions 


c. = c v + e c. 

A ( f -e L ) ( y H - “>) 

H 


(a) 


(5-3) 


s A= s (y H -6 L ) + 6 s (y H -^ } 


+ e s, 


(b) 


are introduced, it can be shown that 


% = | t(3 A°4' H +0iS ^H 6 eH )1+l '' Sfi3 ^H + ^ C ^H' :e H )J ' °*V ’ <S) 


e„ x k* =l[(-s A s l c , +c . c . c., c , - cs Q s Q 1 )i. 
v - a <p H q* a y H e H 0» A Q h Q - 


- ( V<6 h V + c aV H c 0 H °s' ,i 


H H 


(5*4) 


-(s.c, s^, + o A s , c r, s rtl )kl 

A f H 6 ' A e H e “ 


(b) 


o iv = «v'i'=iU t V H V -V*» H « e .«e k - 3 s h 0 s' 0 a) (o) 


If the approximate formula ( I 4 . . 17 ) and the relations (£.1) 

through (5»4) are substituted into the torque equation (3.20(A) 

/ 

it follows that 


% = — ARx i + N 1R „ 1 + N sb „ k 


ARy -*■ T ARz 


(5*5) 
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, a a 3 * 2 
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(b) ( 5 . 6 ) 
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When the torque expression (5«5) and relations (4»5) are 
substituted into equations ( £4. . Z4.) 1 , equations (4*3) give the, 
variational equations for the perturbed problem. 

As in the analysis presented in [M.R. ,19731, it is as- 
sumed that the precession rate XI and the aerodynamic torque 
-are of the same order of magnitude. Let <c represent the 
ratio of the magnitude of aerodynamic torque to the rotational 
energy of the orbiting body about its center of mass and define 

*1 ' ( x ^ ; y k ) = 1 X 1 ( x^ y k ) , 1=1, 2, ...» 6 (5 Sf)'- 

G 1 

where j=l,2 , ..., 6, k=l,2, and y^ y 2 represent the fast varia- 

V, 

bias M, v, respectively. Then the variational equations (4*3) 
take the form 

*\ 

^ l^i ^ ^ * i = l , 2 , • « * , 6 (5*8) 

It is to be noted that the perturbing functions X^ ' (x^; y^) ,i=l, 
2 ,..., 6, are periodic functions of each y k ,k=l,2, with period 
2 tc* Thus, the new dynamical system given by (5*8) 
and, (see [A.R. , 1971, Eqs .( 3 *3(b) ) , ( 3 . 8 ) ] and [F.R., 1972>Eq . 
(2.5( c)) ])* 

y x = M = n ( a) ( 5.9) 

1 All equations designated (4»ij)#i>j> nonnegative integers, 
refer to equations given in [F.R.,1972]. 

When A = B, we have k = 0, q = 0,/3= tx/ 2, v = 0,y= 0, 
y\ o (x/2>0) = 1^(0, k) = 0, _&■ = u, and - 0. 
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J 2 = -u*=4 h = -*=| 

can be treated by the method of averaging, 
A transf ormation 

+ e R^’(x ;y k ), i=l,2, ...,6 


(b) (5.9) 


(5.10) 


to new variables is introduced, so that, for suitable functions 
R^'(x^Jy^.) the fast variables y^, k=l,2> are eliminated from 
the transformed, dynamical system to the first order in ^ 

If relations [F.R.,157l,Eq.{1.13( a) ) ] are used, the transform- 
ed ( the averaged) variational equations take the form 


2x 2 % 


x 1 =^ 1 (X) [ f X' i '(x J iy k )dE d^ H , 1=1,2, ...,6 
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dft 


( j) 
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are introduced for economy of notation. After an integration 
over M and from 0 to 2-76 is perfomed, the averaged differ- 
ential equations (5»H) take the form 


y H * - 1 v+v cot a H ■v [I WL 


c A [k„+ek c, , ]cot Q c , 

2hp4 ° 7 o ( f h -oj) e H S' 
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14 10 9 H 0 H 9 9 
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2 14 10 ( - w) e H 0 H 9 
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2 
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2 (y H -o>) 1 - 00 ) 9 h e H 9 

The variables ( will be used to represent 

the first order secular part of the variables in the remaind- 

er of this section 
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Equations (5»13) represent the complete set of equations 
for determining the first order secular rotational motions of 
a uniaxial body under the influence of the restoring aerodyna- 
mic torque. The precession of the orbital plane and the re- 
gression of perigee are both considered sihee Si and lo may 

o 

have nonzero values. The above equations suggest that there 
are long-term secular changes in the angular momentum vector 
h as well as in the rest 6f the variables. The integrals 
k ,m=0, . . . ,ll|., must be evaluated separately and they are func- 

4 

tions of both __ and oO and the atmospheric density p . 

' H / a 

In the following paragraph, the special case of a uniaxial 

body moving along a circular orbit will be examined. 

In case of a circular orbit, the eccentricity e is zero. 

If we further assume that the molecular density of air p • is 

a constant, it is found that the nonzero k ' s are 

m 

k 0 = f a’ k 3 = k 5 = **8 = *10 =|/°a 


and the associated equations of motion are 
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(d) 

(e) (£.1$) 

h = 0 (f) 

where V Q = y^/p* . 

Equations (£. 15 ) show that there is no secular variation 
in h and 0' to the first order, and if comparison is made be- 
tween equations ( 5 * 15 ) and the associated averaged equations 
for perturbation due togravity-gradient torque given in [2], 
it can be seen that they differ only by a constant coefficient. 
In the case of ideally inelastic collisions between the air 
molecules and the surface of the satellite or if the satellite 
is also symmetric about a plane which is perpendicular to the 
axis of symmetry, we have = 0. Then it may be concluded 
that, in this special case, there is no secular change caused 
by the restoring aerodynamical torque. 

6 . First-Order, Averaged Differential Equations for a Uniaxial 
Body with a Surface of Revolution Subjected to Dissipative 
Aerodynamic Torque 

The dissipative aerodynamic torque and an approximation 
to it have been derived in Sections 3 and ij. of this report. 
Suppose that the body possesses a surface of revolution with 
respect to the a' -axis and take the semi-constrained system 
as the reference system. We can write 
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-AD 


^ C f a V 0 f-C n ^n + (-C U 


^y ,,+c 2 3 s £ 


oO 


,) 


I > I 


+(c 32 3 S V ° J y‘ 1 -C 33 » } — 1 1 3 

If the equations of transformation 
(6.1) becomes 

N = 1 co V [(-C n1 ^ - q s 

-AD 2 ' a 0 11 x 23 


(2.1) are used. 



( 6 . 1 ) 

expression 


+(- C H^ yl + C 23 s ^ ^ cOg | ) ^ 


( 6 . 2 ) 


^- C 32 S Sv s ^x- 


The equations of transformation [A.R. , 197l,Eq .( 4*3( a) ) ] 
can be used to transform (6.2) from the body-fixed system to 
the angular momentum system. The trigonometric functions of^ 
and £ can be expressed in terms of the angles yj ,9 ,9 , d> , 
6‘ , cp' through equations (5.4(a)) and identities 

Sr c A * i* • ( e v x k* ) (a) 

0 V $ y V 

(6.3) 

Sr S k* = i* X (e X k') (b) 

f V r V v 

If then relations [A.R. ,197l,Eqs . ( 4«3( a) ) , ( 4*4( a) ) 3 * (5*3) and 
( 5*4) are used, it is found that 
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s r c = 1 fs . ( -c , a , , - s. c, , c Ql ) 

^V/^Vf A 8 

+C A (c ^ h C 0 h C C Q' " S 0 h C ^ ,S 6« - S ^ H C e H % ,)] 

s r s „ = l[s.(-c , c ,,+ s , s , c « ) (6 

f V Z 3 V f L ‘ft H ^ ft 6 

+ c A ( " c d3 c 0 s ^‘ c a> +s fi 3 j* S a. " s , C Q 3 , ,)] ( b > 

A °H ^ 8 0 H ft 0 ft H 0 H ^ 

The components of h in the body-fixed system are 


h x' 

= h £ 3 f 1 s e' 

( a) 
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= h c t' 3 o’ 

(b) 

h z« 

= h °e> 
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It follows that, for a uniaxial body. 
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( a) 


oJ . = b 

y' a 

C ^ ,S Q« 

(b) 

(6.6) 

cJ , = b 
z ' C 

C 0' 

( c) 


since h x , = A <*J > x , , h 

V 

= A - 3 In , — C C'O * • 

y* z 1 z 1 




Through the combined use of fA.R. , 1571, Eq . ( 4»3( a) ) ] > (6*4-) 
and (6.6), equation (6.2) can now be expressed in the form 
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N._ = N 1 + N Ari i + N at , k 
-AD —AD xr ADy •* AD 2 - 


(6.7) 
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It is to be noted that the dissipative torque is a 

continuous function of the slow variables i=l# . . . , 6, and 

the fast variables S and and it is also a periodic function 

1 H 
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of each of f5 and 4 with a period of 2n . If the torque ex- 

H 

pression (6.7) and relations ( 5*9) are substituted into 
[P .R. »1972, Eqs . ( 4»4) [P t 1972 s Eqs . ( 4»3) ] gives the varia- 
tional equations for the perturbed problem. Again* the equa- 
tions, of transformation (5«10) are introduced and the variables 
^ H ,e ' ’t l,!l are use< ^ t0 represent their first order 
secular parts. The associated equations for the secular motions 
(5»H) under the influence of dissipative aerodynamic torque 
become 

Vs = -'V +a 9 • co Vv H ) -k 
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' e ' = - £11) « e , c 9 , 
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It may be seen from equations (6.9) that all six variables 
have secular, long-term variations due to the presence of the 
dissipative aerodynamic torque, N • The terms which involve 
k are more important than those containing k Q and k n since 

0 0 7 

kg and k^ are smaller quantities than k^. 

The problem associated with a circular orbit and constant 
air density can be obtained readily by letting e = 0 and ^ 

= k = O. It — s> . Thft rl -1 •Pf’ATArvhl 

8 9 

equations for this case are 


a 


constant, so that k„ = k„ = 0, k = p . The differential 

0 I a 
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* s 9° s v . 
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Equations (6.10(d)) and (6.10(f)) can be integrated di- 
rectly. It is found* from (6.10), that 


N c t 

tan Ql = tan Ql e 


( 6 . 11 ) 


and from (6.10(f)), that 
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( 6 . 12 ) 


where 


N 0 = °l°a V o ( C 33 - £ll) 


(6.13) 



k° 


and 9' and h are the initial values of 9' and h, respectively, 
o o 

Prom equations (6.11) and (6.12) * it can be seen that* 
under the influence of dissipative aerodynamic torque* the 
magnitude of the angular momentum vector decreases exponential- 
ly, approaching zera in the limit. However* the angle between 
the x* -axis and the angular momentum vector may either de- 
crease or increase with time* approaching either zero or x/2 
as a limiting value* respectively* depending on whether N 0 is 
negative or positive. If and are of the same order of 
magnitude* equation (6.11) indicates that the body will event- 
ually spin about the axis of the maximum moment of inertia. 

In case 9' is zero, 9* will remain zero and hence 9' is a 
o 

constant of motion. 
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